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Abstract. We enumerate all totally real number fields F with root 
discriminant Sf < 14. There are 1229 such fields, each with degree 
[F : Q] < 9. 

In this article, we consider the following problem. 

Problem 1. Given B £ R>o, enumerate the set NF(B) of totally real number 
fields F with root discriminant Sf < B, up to isomorphism. 

To solve Problem 1, for each n 6 Z>o we enumerate the set 

NF(n, B) = {F e NF(B) : [F : Q] = n} 

which is finite (a result originally due to Minkowski). If F is a totally real field 
of degree n — [F : Q], then by the Odlyzko bounds [57], we have Sf > 47re 1+7 — 
<3(n~ 2 / 3 ) where 7 is Euler's constant; thus for B < 4ne 1+1 < 60.840, we have 
NF(n, B) = for n sufficiently large and so the set NF(B) is finite. Assuming 
the generalized Riemann hypothesis (GRH), we have the improvement Sp > 
87re 7+Tr / 2 — 0(log~ 2 n) and hence NF(B) is conjecturally finite for all B < 
87re 7+7r/2 < 215.333. On the other hand, for B sufficiently large, the set NF(B) 
is infinite: Martin [53] has constructed an infinite tower of totally real fields 
with root discriminant i5 « 913.493 (a long-standing previous record was held 
by Martinet [25] with S w 1058.56). The value 

liminf mm{S F : F € NF(n, B)} 

n — 'oc 

is presently unknown. If B is such that #NF(B) = 00, then to solve Problem 1 
we enumerate the set NF(B) = \J n NF(n, B) by increasing degree. 

Our restriction to the case of totally real fields is not necessary: one may 
place alternative constraints on the signature of the fields F under considera- 
tion (or even analogous p-adic conditions). However, we believe that Problem 1 
remains one of particular interest. First of all, it is a natural boundary case: by 
comparison, Hajir-Maire [14115] have constructed an unramified tower of totally 
complex number fields with root discriminant w 82.100, which comes within a 



factor 2 of the GRH-conditional Odlyzko bound of 87re 7 « 44.763. Secondly, 
in studying certain problems in arithmetic geometry and number theory — for 
example, in the enumeration of arithmetic Fuchsian groups [21j and the compu- 
tational investigation of the Stark conjecture and its generalizations — provably 
complete and extensive tables of totally real fields arc useful, if not outright 
essential. Indeed, existing strategies for finding towers with small root discrim- 
inant as above often start by finding a good candidate base field selected from 
existing tables. 

The main result of this note is the following theorem, which solves Problem 
1 for 6 = 14. 

Theorem 2. We have #NF(U) = 1229. 

The complete list of these fields is available online [35]; the octic and nonic 
fields (n = 8,9) are recorded in Tables 3-4 in §4, and there are no dectic fields 
(iV.F(14, 10) = 0). For a comparison of this theorem with existing results, see 
§1.2. 

The note is organized as follows. In §1, we set up the notation and back- 
ground. In §2, we describe the computation of primitive fields F € NF(14); we 
compare well-known methods and provide some improvements. In §3, we discuss 
the extension of these ideas to imprimitive fields, and we report timing details 
on the computation. Finally, in §4 we tabulate the fields F. 

The author wishes to thank: Jiirgen Kliiners, Noam Elkies, Claus Fieker, Ki- 
ran Kedlaya, Gunter Malle, and David Dummit for useful discussions; William 
Stein, Robert Bradshaw, Craig Citro, Yi Qiang, and the rest of the Sage devel- 
opment team for computational support (NSF Grant No. 0555776); and Larry 
Page and Helen Read for their technical assistance. 

1 Background 
1.1 Initial bounds 

Let F denote a totally real field of degree n = [F : Q] with discriminant dp and 
root discriminant 5p = d% n . By the unconditional Odlyzko bounds [27] (see also 
Martinet [24]), if n > 11 then S F > 14.083, thus if F G iVF(14) then n < 10. 



Table 1. Degree and Root Discriminant Bounds 



n 


2 345678 9 10 


Bo 
Bo (GRH) 


> 2.223 3.610 5.067 6.523 7.941 9.301 10.596 11.823 12.985 

> 2.227 3.633 5.127 6.644 8.148 9.617 11.042 12.418 13.736 


A 


30 25 20 17 16 15.5 15 14.5 14 



The lower bounds for 5f in the remaining degrees are summarized in Table 
1: for each degree 2 < n < 10, we list the unconditional Odlyzko bound Bo = 



Bo {n) , the GRH-conditional Odlyzko bound (for comparison only, as computed 
by Cohen-Diaz y Diaz-Olivier [7]), and the bound 5f < A that we employ. 

1.2 Previous work 

There has been an extensive amount of work done on the problem of enumerating 
number fields — we refer to [18] for a discussion and bibliography. 

1. The KASH and PARI groups [16] have computed tables of number fields of 
all signatures with degrees < 7: in degrees 6, 7, they enumerate totally real 
fields up to discriminants 10 7 , 15 • 10 7 , respectively (corresponding to root 
discriminants 14.67, 14.71, respectively). 

2. Malle [53] has computed all totally real primitive number fields of discrim- 
inant dp < 10 9 (giving root discriminants 31.6, 19.3, 13.3, 10 for degrees 
6,7,8,9). This was reported to take several years of CPU-time on a SUN 
workstation. 

3. The database by Kliiners-Malle [T7] contains polynomials for all transitive 
groups up to degree 15 (including possible combinations of signature and 
Galois group); up to degree 7, the fields with minimal (absolute) discriminant 
with given Galois group and signature have been included. 

4. Roblot [30] constructs abelian extensions of totally real fields of degrees 4 to 
48 (following Cohen-Diaz y Diaz-Olivier [B]) with small root discriminant. 

The first two of these allow us only to determine NF (10) (if we also separately 
compute the imprimitive fields); the latter two, though very valuable for certain 
applications, are in a different spirit than our approach. Therefore our theorem 
substantially extends the complete list of fields in degrees 7-9. 

2 Enumeration of totally real fields 
2.1 General methods 

The general method for enumerating number fields is well-known (see Cohen [3] 
§9.3]). We define the Minkowski norm on a number field F by T^a) = \ a i\ 2 
for a G F, where ax, ol%, ■ ■ . ,a n are the conjugates of a in C. The norm T2 gives 
Zp the structure of a lattice of rank n. In this lattice, the element 1 is a shortest 
vector, and an application of the geometry of numbers to the quotient lattice 
Zf/Z yields the following result. 

Lemma 3 (Hunter). There exists a E Zp \ Z such that < Tr(a) < n/2 and 




where 7„_i is the (n— l)th Hermite constant. 



Remark 4- The values of the Hermite constant are known for n < 8 (given by the 
lattices A x , A 2 , A 3 , D 4 , D 5 , E 6 , E 7 , E 8 ): we have 7™ = 1,4/3,2,4,8,64/3,64,256 
(see Conway and Sloane [S]) for n = 1, . . . , 8; the best known upper bounds for 
n = 9, 10 are given by Cohn and Elkies [5]. 

Therefore, if we want to enumerate all number fields F of degree n with 
\dp\ B, an application of Lemma [3] yields a € 7Lp\7L such that 12(a) < C for 
some C € R>o depending only on n, B. We thus obtain bounds on the power 
sums 

n n 

\S k (a)\ =J2 a i ^ T *(«) = E ^ " Cfc/2 ' 
and hence bounds on the coefficients a% € Z of the characteristic polynomial 

i=l 

of a by Newton's relations: 

fc-i 



a 



+ ^ a n-%Sk-i + ka„- k = 0. (1) 



This then yields a finite set NS(n, B) of polynomials f(x) € Z[x] such that every 
F is represented as Q[x]/(f(x)) for some f(x) £ NS(n, B), and in principle each 
f(x) can then be checked individually. We note that it is possible that a as 
given by Hunter's theorem may only generate a subficld Q C Q(a) C F if F is 
imprimitive: for a treatment of this case, see §3. 

The size of the set NS[n,B) is 0(£ n ("+ 2 )/ 4 ) (see Cohen [1 §9.4]), and the 
exponential factor in n makes this direct method impractical for large n or B. 
Note, however, that it is sharp for n — 2: we have NF(2,B) ~ (6/ir 2 )B 2 (as 
i? — > 00), and indeed, in this case one can reduce to simply listing squarefree 
integers. For other small values of n, better algorithms are known: following 
Davenport- Hcilbronn, Belabas [5] has given an algorithm for cubic fields; Cohen- 
Diaz y Diaz-Olivier [7] use Kummer theory for quartic fields; and by work of 
Bhargava [3] , in principle one should similarly be able to treat the case of quintic 
fields. No known method improves on this asymptotic complexity for general n, 
though some possible progress has been made by Ellenberg-Venkatesh [12] . 



2.2 Improved methods for totally real fields 

We now restrict to the case that F is totally real. Several methods can then be 
employed to improve the bounds given above — although we only improve on the 
implied constant in the size of the set NS(n, B) of examined polynomials, these 
improvements are essential for practical computations. 



Basic bounds. From LemmaO we have < a„_i = — Tr(a) < \n/2\ and 



1 -2 Iw^s 1 ^ !\ „2 7»-l (B 



a n -2 = -a„_ : - -T 2 (a) > - 1 - - a n _ Y 



l/(n-l) 



For an upper bound on a„_ 2 , we apply the following result. 
Lemma 5 (Smyth |32j). If ^ is a totally positive algebraic integer, then 

Tr( 7 ) > 1.7719[Q( 7 ) : Q] 
unless j is a root of one of the following polynomials: 

x-l,x 2 -3x+l,x 3 -5x 2 + 6x-l,x 4 -7x 3 + 13x 2 -7x+l,x 4 ~7x 3 + Ux 2 -8x + l. 

Remark 6. The best known bound of the above sort is due to Aguirre-Bilbao- 
Peral [J, who give Tr(7) > 1.780022[Q(7) : Q] with 14 possible explicit excep- 
tions. For our purposes (and for simplicity), the result of Smyth will suffice. 

Excluding these finitely many cases, we apply Lemma[5]to the totally positive 
algebraic integer a 2 , using the fact that T 2 (a) = Tr(a 2 ), to obtain the upper 
bound a„_ 2 < <4_i/2 - 0.88595n. 



Rolle's theorem. Now, given values a n _i, a rl _ 2 , • ■ ■ , cin-k f° r the coefficients of 
f{x) for some k > 2, we deduce bounds for a n -k-i using Rolle's theorem — this 
elementary idea can already be found in Takeuchi [33] and Kliiners-Malle [HI 
§3.1]. Let 

f(n-i)r x \ 
fi( x ) = — f TT~ = 9i( x ) + a n-i 

[n — 1)\ 

for i = 0, . . . , n. Consider first the case k = 2. Then 

, , n(n - l)(n - 2) , (n-l)(n-2) , . 
g 3 (x) = -i £ J -x 3 + ± ^ '-dn-xx 2 + (n - 2)a„_ 2 a;. 

Let /?i < /?2 denote the roots of / 2 (x). Then by Rolle's theorem, 

/3(A) = 03 (A) + «n- 3 > and f 3 (f3 2 ) = 03(^2) + a„_ 3 < 

hence — 53 < a„_3 < — gz{pi). In a similar way, if /3j < • • • < denote 
the roots of fk(x), then we find that 

- mm gk+i(Pi } ) < On-k-i < ~ max gk+i(j3f '). 

l<t<k l<i<k 
i^k (2) i=k (2) 



Lagrange multipliers. We can obtain further bounds as follows. We note that 

(k) / (k) 

if the roots of / are bounded below by /3q (resp. bounded above by /3j, +1 ), then 

fM k) ) =g k {P { a y ) + a n -k >0 

(with a similar inequality for p^)-^, and these combine with the above to yield 
~ nJ^jV, , 9k+i(Pi) < a n -k-i < - max g k+1 {l3f ] ). (2) 

i^k (2) i=k (2) 

(k) (k) 

We can compute /?q , by the method of Lagrange multipliers, which 
were first introduced in this general context by Pohst [29] (see Remark [7]). The 
values a„_x, ■ ■ ■ , a n -k G ^ determine the power sums Sj for i = 1, . . . , k by 
Newton's relations |T]). Now the set of all a; = (a;,) E R" such that Si(a;) = Sj is 
closed and bounded, and therefore by symmetry the minimum (resp. maximum) 
value of the function x n on this set yields the bound /3q (resp. P^li)- By the 
method of Lagrange multipliers, we find easily that if a; € W 1 yields such an 
extremum, then there are at most k — 1 distinct values among x\, . . . , x n -i, 
from which we obtain a finite set of possibilities for the extremum x. 

For example, in the case k = 2, the extrema are obtained from the equations 

(n - l)xi + x n = s 1 = — a„_i and (n - l)xf + x 2 n — s 2 = a\_ 1 - 2a„_ 2 
which yields simply 



P ( o 2 \p { 3 2) = -\ ± i» - : ; 2(^1 ■ — T j a?l _ 2 j . 

(It is easy to show that this always improves upon the trivial bounds used by 
Takeuchi [33 .) For k — 3, for each partition of n — 1 into 2 parts, one obtains 
a system of equations which via elimination theory yield a (somewhat lengthy 
but explicitly given) degree 6 equation for x n . For k > 4, we can continue 
in a similar way but we instead solve the system numerically, e.g., using the 
method of homotopy continuation as implemented by the package PHCpack 
developed by Verschelde [34] : in practice, we do not significantly improve on these 
bounds whenever k > 5, and even for k = 5, if n is small then it often is more 
expensive to compute the improved bounds than to simply set (3^ = ^ 

Remark 7. Pohst's original use of Lagrange multipliers, which applies to number 
fields of arbitrary signature, instead sought the extrema of the power sum Sk+i to 
bound the coefficient a n -k-i- The bounds given by Rolle's theorem for totally 
real fields are not only easier to compute (especially in higher degree) but in 
most cases turn out to be strictly stronger. We similarly find that many other 
bounds typically employed in this situation (e.g., those arising from the positive 
definiteness of T 2 on Z[a\) are also always weaker. 
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2.3 Algorithmic details 



Our algorithm to solve Problem 1 then runs as follows. We first apply the basic 
bounds from §2.2 to specify finitely many values of a n _i,a n _2- For each such 
pair, we use Rolle's theorem and the method of Lagrange multipliers to bound 
each of the coefficients inductively. Note that if k > 3 is odd and a„_i = a„_3 = 
• • • = a n _(fc_2) = 0, then replacing a by —a we may assume that a n -u > 0. 

For each polynomial / £ NS(n, B) that emerges from these bounds, we test 
it to see if it corresponds to a field F £ NF(n, B). We treat each of these latter 
two tasks in turn. 

Calculation of real roots. In the computation of the bounds we use 
Newton's method to iteratively compute approximations to the roots /3. , using 
the fact that the roots of a polynomial are interlaced with those of its derivative, 
i.e. /3^7 X) < /?f 3 < 0i k _1) for i = l,...,k. Note that by Rolle's theorem, we 
will either find a simple root in this open interval or we will converge to one 
of the endpoints, say = 0+ ^, and then necessarily (3^ = ft^ 1 ' as well, 
which implies that fk(x) is not squarefree and hence the entire coefficient range 
may be discarded immediately. It is therefore possible to very quickly compute 
an approximate root which differs from the actual root /3 J - fc+1 ^ by at most some 
fixed e > 0. We choose e small enough to give a reasonable approximation but 
not so small as to waste time in Newton's method (say, e = 1CP 4 ). We deal with 
the possibility of precision loss by bounding the value gk+i{!3^) in @ using 
elementary calculus; we leave the details to the reader. 

Testing polynomials. For each / <E NS(n,B), we test each of the following 
in turn. 

1. We first employ an "easy irreducibility test": We rule out polynomials / 
divisible by any of the factors: x, x ± 1, x ± 2, x 2 ± x — 1, x 2 — 2. In the 
latter three cases, we first evaluate the polynomial at an approximation to 
the values (1 ± v5)/2, V%, respectively, and then evaluate / at these roots 
using exact arithmetic. (Some benefit is gained by hard coding this latter 
evaluation.) 

2. We then compute the discriminant d = disc(/). If d < 0, then / is not a real 
separable polynomial, so we discard /. 

3. If F = Q[a] = Q[x]/(f(x)) £ NF(n,B), then for some a £ Z we have 
Bo{n) n < dp = d/a 2 < B n where Bo is the Odlyzko bound (see §1). 
Therefore using trial division we can quickly determine if there exists such 
an a 2 | d; if not, then we discard /. 

4. Next, we check if / is irreducible, and discard / otherwise. 

5. By the preceding two steps, an a-maximal order containing Z[a] is in fact 
the maximal order Zf of the field F. If disc(ZF) = c?f > B, we discard /. 

6. Apply the POLRED algorithm of Cohen-Diaz y Diaz 0: embed Zf C 1" 
by Minkowski (as in §1.1) and use LLL- reduction [20] to compute a small 



element a re d € If such that Q(a) — Q(a rc d) = F. Add the minimal polyno- 
mial f le d(z) of a re d to the list NF(n,B) (along with the discriminant (If), 
if it does not already appear. 

We expect that almost all isomorphic fields will be identified in Step 6 by 
computing a reduced polynomial. For reasons of efficiency, we wait until the 
space NS(n, B) has been exhausted to do a final comparison with each pair of 
polynomials with the same discriminant to see if they are isomorphic. Finally, 
we add the exceptional fields coming from Lemma [5] if relevant. 

Remark 8. Although Step 1 is seemingly trivial, it rules out a surprisingly sig- 
nificant number of polynomials / — indeed, nearly all reducible polynomials are 
discarded by this step in higher degrees. Indeed, if T 2 (f) = J2i a i (where a« are 
the roots of /) is small compared to deg(/) = n, then / is likely to be reducible 
and moreover divisible by a polynomial g with T 2 {g) also small. It would be 
interesting to give a precise statement which explains this phenomenon. 

2.4 Implementation details 

For the implementation of our algorithm, we use the computer algebra system 
Sage [31], which utilizes PARI [28] for Steps 4-6 above. Since speed was of the 
absolute essence, we found that the use of Cython (developed by Stein and Brad- 
shaw) allowed us to develop a carefully optimized and low-level implementation 
of the bounds coming from Rollc's theorem and Lagrange multiplier method 
O We used the DSage package (due to Qiang) which allowed for the distribu- 
tion of the compution to many machines; as a result, our computational time 
comes from a variety of processors (Opteron 1.8GHz, Athlon Dual Core 2.0GHz, 
and Celeron 2.53GHz), including a cluster of 30 machines at the University of 
Vermont. 

In low and intermediate degrees, where we expect comparatively many fields, 
we find that the running time is dominated by the computation of the maximal 
order (Step 5), followed by the check for irreducibility (Step 4); this explains the 
ordering of the steps as above. By contrast, in higher degrees, where we expect 
few fields but must search in an exponentially large space, most of the time is 
spent in the calculation of real roots and in Step 1. Further timing details can 
be found in Table 2 in §3.2. 

3 Imprimitive fields 

In this section, we extend the ideas of the previous section to imprimitive fields 
F, i.e. those fields F containing a nontrivial subfield. Suppose that F is an 
extension of E with [F : E] = m and [E : Q] = d. Since S F > S E , if F e NF(B) 
then E G NF(B) as well, and thus we proceed by induction on E. For each such 
subfield E, we proceed in an analogous fashion. We let 

f(x) = x m + arn-ix" 1 ^ 1 H h a x x + a 

be the minimal polynomial of an element agZf with F = E(a) and ai € 1i E . 



3.1 Extension of bounds 



Basic bounds We begin with a relative version of Hunter's theorem. We denote 
by -Boo the set of infinite places of E. 

Lemma 9 (Martinet [26] ). There exists a £ 7Lp\%E such that 

U«) < I £ |„ «)| 2 +7-J ( J|L) " M . (3) 

The inequality of Lemma[H]remains true for any element of the set /iecx + Ze, 
where he denotes the roots of unity in E. This allows us to choose TTp/E a = 
— a m _i among any choice of representatives from Z^/mZE (up to a root of 
unity); we choose the value of a m _i which minimizes 

^ |o-(Tr F /j5a)| = ^2 v(am-i) 2 , 

which is a positive definite quadratic form on Z^; such a value can be found 
easily using the LLL-algorithm. 

Now suppose that F is totally real. Then ^2 a£E \ <7 ( a m-i)\ 2 — ^E/Q^m-ii 
and we have 2 d or \m d /2] possibilities for a m _i, according as m — 2 or otherwise. 
For each value of a m _i, we have T2{a) £ Z bounded from above by Lemma[5] 
and from below by Lemma O since Ti e /q(c( 2 ) = T 2 (a) > 1.7719n. If we denote 
Trpy^a 2 = t 2 , then by Newton's relations, we have £2 = a m-i — 2a m _2, and 
hence Tt^/q^ = ?2(a) and £2 = a m-i ( m °d 2). In particular, t-i £ Ze is totally 
positive and has bounded trace, leaving only finitely many possibilities: indeed, 
if we embed Ze K d by Minkowski, these inequalities define a parallelopiped 
in the positive orthant. 



Lattice points in boxes. One option to enumerate the possible values of ti 
is to enumerate all lattice points in a sphere of radius given by ^ using the 
Fincke-Pohst algorithm [13] , However, one ends up enumerating far more than 
what one needs in this fashion, and so we look to do better. The problem we 
need to solve is the following. 

Problem 10. Given a lattice L c M d of rank d and a convex polytope P of finite 
volume, enumerate the set POL. 

Here we must allow the lattice L to be represented numerically; to avoid 
issues of precision loss, one supposes without loss of generality that dPC\L = 0. 

There exists a vast literature on the classical problem of the enumeration of 
integer lattice points in rational convex polytopes (see e.g., De Loera [TD]), as well 
as several implementations jllll9j . (In many cases, these authors are concerned 
primarily with simply counting the number of lattice points, but their methods 
equally allow their enumeration.) 



In order to take advantage of these methods to solve Problem[Tni we compute 
an LLL- reduced basis 7 = 71 , . . . , 7^ of L, and we perform the change of variables 
4> : M. d — > M. d which maps ji 1— > where is the ith coordinate vector. The 
image 4>(P) is again a convex polytope. We then compute a rational polytope Q 
(i.e. a polytope with integer vertices) containing 4>{P) by rounding the vertices 
to the nearest integer point as follows. For each pair of vertices v,w € P such 
that the line £(v, w) containing v and w is not contained in a proper face of P, 
we round the ith coordinates 4>(v)i down and 4>{w)i up if <fr(v)i < 4>{w)i, and 
otherwise round in the opposite directions. The convex hull Q of these rounded 
vertices clearly contains <j)(P H L), and is therefore amenable to enumeration 
using the methods above. 

We note that in the case where P is a parallelopiped, for each vertex v there 
is a unique opposite vertex w such that the line l(y, w) is not contained in a 
proper face, so the convex hull Q will also form a parallelopiped. 



Coefficient bounds and testing polynomials. The bounds in §2 apply mu- 
tatis mutandis to the relative situation. For example, given a m _i, . . . , s m _t for 
k > 2, for each v S E^,, if we let v(g) denote the polynomial Y^i v Q 3 i) xl f° r 
g{x) = J^i biX 1 G E[x], we obtain the inequality 

~ n J°°}? , v (9k+i)(P- k y) < w(a m _ fc _i) < - max u(£? fc+ i)(/3|^); 

0<«<fc+l 0<2<fc + l 

i^/c (2) t=fc (2) 

here, /3^, . . . denote the roots of v(f k (x)), and /3o^>/^[+i, v are computed in 
an analogous way using Lagrange multipliers. In this situation, we have a rn -k-i 
contained in an honest rectangular box, and the results of the previous subsection 
apply directly. 

For each polynomial which satisfies these bounds, we perform similar tests 
to discard polynomials as in §2.3. One has the option of working always relative 
to the ground field or immediately computing the corresponding absolute field; 
in practice, for the small base fields under consideration, these approaches seem 
to be comparable, with a slight advantage to working with the absolute field. 



3.2 Conclusion and timing 

Putting together the primitive and imprimitive fields computed in §§2-3, we 
have proven Theorem 2. In Table 2, we list some timing details arising from the 
computation. Note that in high degrees (presumably because we enumerate an 
exponentially large space) we recover all imprimitive fields already during the 
search for primitive fields. 



Table 2. Timing data 



n 


2 3 


4 


5 


6 


7 


8 


9 


10 


A(n) 


30 25 


20 


17 


16 


15.5 


15 


14.5 


14 


f 


443 4922 


57721 


244600 3242209 1.7 x 10 r 


1.2 x 10 8 


9.5 x 10 s 


2.2 x 10 9 


Irred / 


418 2523 


27234 


157613 2710965 1.6 x 10 7 


1.1 x 10 s 


9.0 x 10 s 


2.2 x 10 9 


f,d F <B 


418 1573 


5665 


4497 


1288 


4839 


3016 


506 





F 


273 630 


1273 


674 


802 


301 


164 


15 





Total time 


0.2s 2.2s 


26.8s 


lm25s 


17m3s 


2h59m 


ld4.5h 


17d21h 


173d 


Imprim / 





7059 





62532 





239404 


15658 


945866 


Imprim F 





702 





420 





100 


6 





Time 




4m22s 




8m38s 




lh56m 


16m53s 


llh27m 


Total fields 


273 630 


1578 


674 


827 


301 


164 


15 






4 Tables of totally real fields 

In Table 3, we count the number of totally real fields F with root discriminant 
Sp < 14 by degree, and separate out the primitive and imprimitive fields. We also 
list the minimal discriminant and root discriminant for n < 9. The polynomial 

x 10 - 11a; 8 - 3a; 7 + 37a; 6 + 14a; 5 - 48a; 4 - 22a; 3 + 20a; 2 + 12a; + 1 

with d F = 443952558373 = 61 2 397 2 757 and S F » 14.613 is the dectic totally real 
field with smallest discriminant that we found — the corresponding number field 
(though not this polynomial) already appears in the tables of Kluners-Malle 
p"7] and is a quadratic extension of the second smallest real quintic field, of 
discriminant 24217. It is reasonable to conjecture that this is indeed the smallest 
such field. 



Table 3. Totally real fields F with 6 F < 14 



n = [F : Q] 


#NF(n, 14) 


Primitive F 


Imprimitive F 


Minimal dp 


Minimal 8f 


2 


59 


59 





5 


2.236 


3 


86 


86 





49 


3.659 


4 


277 


117 


160 


725 


5.189 


5 


170 


170 





14641 


6.809 


6 


263 


104 


159 


300125 


8.182 


7 


301 


301 





20134393 


11.051 


8 


62 


19 


43 


282300416 


11.385 


9 


11 


6 


5 


9685993193 


12.869 


10 











443952558373? 


14.613? 


Total 


1229 


862 


367 





In Tables 4-5, we list the octic and nonic fields F with 5p < 14. For each 
field, we specify a maximal subfield E by its discriminant and degree — when more 
than one such subfield exists, we choose the one with smallest discriminant. 



Table 4. Octic totally real fields F with S F < 14 



dp 






/ 


[E : Q 


!] d E 


282300416 




x" - Ax' + 14x* - 8a:* - 12x 3 + 7a 2 + 2a - 1 


A 


2624 


309593125 


x 8 


- 4x 7 


- a 6 + 17a 5 - 5a; 4 - 23a; 3 + 6a 2 + 9a - 1 


A 


725 


324000000 






a; 8 - 7a 6 + 14a; 4 - 8a; 2 + 1 


A 


1125 


410338673 


x 8 


-x 7 - 


7a 6 + 6x 5 + 15a 4 - 10a 3 - 10a; 2 + 4x + 1 


A 


4913 


432640000 


x 8 


- 2x 7 


- 7x 6 + 16x 5 + 4x 4 - 18x 3 + 2a 2 + 4a - 1 


A 


1600 


442050625 


x 8 - 


2x 7 - 


12a 6 + 26a 5 + 17x 4 - 36x 3 - 5x 2 + 11a - 1 


4 


725 


456768125 


X s - 


- 2x 7 - 


- 7a 6 + 11a 5 + 14x 4 - 18x 3 - 8x 2 + 9x - 1 


4 


725 


483345053 




x 8 - x 7 - 7x 6 + 4a; 5 + 15x 4 - 3a; 3 - 9a; 2 + 1 


1 


1 


494613125 


X 


8 -x 7 


- 7a 6 + 4a 5 + 13a 4 - 4x 3 - 7a 2 + a + 1 


4 


725 


582918125 


x 8 


-2x 7 


- 6a 6 + 9a 5 + 11a 4 - 9x 3 - 6a 2 + 2x + 1 


4 


725 


656505625 


x 8 


- 3a; 7 


- 4x 6 + 13x 5 + 5x 4 - 13x 3 - 4x 2 + 3x + 1 


4 


725 


661518125 


x 8 


-x 7 - 


- 7x 6 + 5x 5 + 15a 4 - 7a 3 - 10a 2 + 2x + 1 


2 


5 


707295133 




x 8 - 8 


a 6 - 2x 5 + 19x 4 + 7x 3 - 13a 2 - Ax + 1 


1 


1 


733968125 


x 8 - 


- 2a; 7 - 


- 6x 6 + 10a 5 + 11a 4 - 11a 3 - 7x 2 + 2a + 1 


2 


5 


740605625 


x 8 


-x 7 - 


9a 6 + 8x 5 + 21a 4 - 12a 3 - 14x 2 + 4x + 1 


4 


725 


803680625 


x 8 - 


2a; 7 - 


9x 6 + 12x 5 + 22a 4 - 24a 3 - 14x 2 + 14a - 1 


4 


725 


852038125 




x 8 - 


■ 10x 6 - 5a 5 + 17a 4 + 5a 3 - 10x 2 + 1 


4 


725 


877268125 


X s 


- 3a; 7 


- 6a 6 + 20x 5 + 5a 4 - 25a 3 - a 2 + 7x + 1 


4 


725 


898293125 


x 8 


-a; 7 - 


- 9x 6 + 10a 5 + 15a 4 - 10a 3 - 9x 2 + x + 1 


4 


725 


1000118125 


x 8 


- 3a; 7 


- 4x 6 + 14x 5 + 5x 4 - 19x 3 - a 2 + 7a - 1 


2 




1024000000 






a 8 - 8x 6 + 19x 4 - 12x 2 + 1 


4 


1600 


1032588125 




x 8 - 9a; 6 - 2a; 5 + 23a; 4 + 9a; 3 - 17a; 2 - 9a; - 1 


2 




1064390625 






a 8 - 13a 6 + 44x 4 - 17a 2 + 1 


4 


725 


1077044573 


X s 


-a; 7 - 


- 8a 6 + 8x 5 + 16a 4 - 17a 3 - 2a 2 + 5a - 1 


1 


1 


1095205625 


x 8 


- 3a; 7 


- 5x 6 + 18x 5 + 2x 4 - 23x 3 + 2a 2 + 8a - 1 


2 


5 


1098290293 


x 8 


- 3a; 7 


- 4a 6 + 16x 5 + x 4 - 23x 3 + 7x 2 + 5a - 1 


1 


1 


1104338125 


x 8 - 


-2a; 7 - 


8a 6 + 15x 5 + 17x 4 - 31x 3 - 9a 2 + 17a - 1 


4 


725 


1114390153 




a; 8 - 


- 8x 6 - 2a 5 + 16x 4 + 3a 3 - 10a 2 + 1 


1 


1 


1121463125 




a; 8 - 3a; 7 - 4a; 6 + 15a; 5 + 2a; 4 - 18a; 3 + 5a; + 1 


2 


5 


1136700613 


X 


8 - a; 7 


- 7a 6 + 4a 5 + 14a 4 - 4x 3 - 8a 2 + a + 1 


1 


1 


1142440000 


X s 


- 3a 7 


- 5x 6 + 15x 5 + 8x 4 - 15x 3 - 5x 2 + 4x + 1 


4 


4225 


1152784549 


X s 


- 4a 7 


- a 6 + 15a 5 - 3x 4 - 16x 3 + 4a 2 + 4a - 1 


4 


1957 


1153988125 


a 8 - 


-2a; 7 - 


- 7a 6 + 11a 5 + 12x 4 - 16x 3 - 5x 2 + 6x - 1 


4 


2525 


1166547493 


a ! 


! - a 7 


- 7x 6 + 6a 5 + 14x 4 - 9a 3 - 9a 2 + 3a + 1 


1 


1 


1183423341 


x 8 


-x 7 - 


8a 6 + 9a 5 + 17x 4 - 20a 3 - 8x 2 + lOx - 1 


4 


1957 


1202043125 




x 8 -3 


a 7 - 4x 6 + 16x 5 - 21x 3 + 9a 2 + 2a - 1 


2 


5 


1225026133 


x 8 


- 3a 7 


- 4x 6 + 18a 5 - 6a 4 - 17a 3 + 9a 2 + 2x - 1 


1 


1 


1243893125 


x 8 -x 7 


- 8x 6 + 3x 5 + 18x 4 - x 3 - 12x 2 - 2a + 1 


2 


5 


1255718125 


x 8 - 


2a 7 - 


8a 6 + 19x 5 + 10a 4 - 41x 3 + 13a 2 + 10a - 1 


4 


725 


1261609229 


x 8 


- 2a 7 


- 6x 6 + 12x 5 + 9x 4 - 19x 3 - a 2 + 6a - 1 


1 


1 


1292203125 


x 8 


- 4a; 7 


- x 6 + 17a 5 - 6x 4 - 21a 3 + 6a 2 + 8x + 1 


4 


1125 


1299600812 


x 8 - 


- 2x 7 - 


- 6x 6 + 10a 5 + 12a 4 - 13a 3 - 8x 2 + 3a + 1 


1 


1 


1317743125 


x 8 


-x 7 - 


8a 6 + 7a 5 + 19x 4 - 14a 3 - 12a 2 + 8x - 1 


2 


5 


1318279381 


X 


8 - a; 7 


- 7a 6 + 5a 5 + 14a 4 - 6x 3 - 9a 2 + a + 1 


1 


1 


1326417388 


x 8 - 


-2a; 7 - 


- 6x 6 + 10a 5 + 12a 4 - 13a 3 - 9x 2 + 4a + 2 


4 


2777 


1348097653 


x 8 - 


- 2a; 7 - 


- 6x 6 + 11a 5 + 11a 4 - 17a 3 - 6x 2 + 6a + 1 


1 


1 


1358954496 






a 8 - 8x 6 + 20x 4 - 16x 2 + 1 


4 


2048 


1359341129 




a; 8 - 


8x 6 - a 5 + 18x 4 + 2x 3 - 12a 2 - a + 2 


1 


1 


1377663125 




a; 8 - 


- 12x 6 + 33x 4 - 5x 3 - 22a 2 + 5x + 1 


4 


725 


1381875749 


x 8 


- 3a; 7 


- 4a 6 + 14x 5 + 4a 4 - 18a 3 + x 2 + 5a - 1 


1 


1 


1391339501 




a; 8 - 3a 7 - 4a; 6 + 15a; 5 + 4a; 4 - 22a 3 + 9a - 1 


1 


1 


1405817381 




a; 8 - 9a 6 - a; 5 + 20a 4 + 6a 3 - 12x 2 - 7a - 1 


1 


1 


1410504129 




a; 8 - 


9a 6 - a 5 + 22x 4 + x 3 - 15x 2 - x + 1 


4 


3981 


1410894053 


x 8 


- 2a 7 


- 6x 6 + 9a 5 + 12x 4 - llx 3 - 8x 2 + 3x + 1 


1 


1 


1413480448 


x 8 


- 4a 7 


- 2x 6 + 16x 5 - x 4 - 16x 3 + 2a 2 + 4a - 1 


4 


2048 


1424875717 


8 7 

x — x 


- 7x 6 + 5a 5 + 15x 4 - 6a 3 - 10x 2 + a + 1 


1 


1 


1442599461 


x 8 


- 3a 7 


- 4x 6 + 15x 5 + 4x 4 - 21x 3 - 2x 2 + 8x + 1 


4 


7053 


1449693125 


a 8 - 


-a 7 - 


9x 6 + 10a 5 + 20a 4 - 20a 3 - 14x 2 + 11a + 1 


2 


5 


1459172469 


x 8 


- 4a 7 


- a 6 + 17a 5 - 6x 4 - 21x 3 + 8a 2 + 6a - 1 


4 


1957 


1460018125 


x 8 - 


- 3a; 7 - 


- 5x 6 + 13a 5 + 11a 4 - 14a 3 - 10a 2 + a + 1 


4 


2525 


1462785589 


x 8 - 


-2a; 7 - 


- 6a 6 + 11a 5 + 10x 4 - 17x 3 - 3x 2 + 6x - 1 


1 


1 


1472275625 


x 8 - 


3a 7 - 


6x 6 + 19x 5 + 13a 4 - 35a 3 - 12x 2 + 13a - 1 


4 


725 



Table 5. Nonic totally real fields F with S F < 14 











/ 








[E:Q 


>] d E 


9685993193 




X' 


1 - 9a; 


7 + 24a; 5 - 2a; 4 - 


20a; 3 


+ 3a; 2 + 5a,- - 1 




1 


1 


11779563529 




x 9 - 


9a; 7 - 


2a; 6 + 22a; 5 + 5a; 


1 - 17a; 3 - 4a; 2 + 4a; + 1 




1 


1 


16240385609 


ar 




5 -9x 


7 + 4a; 6 + 26a; 5 - 


2a; 4 - 


- 25a; 3 - a; 2 + 7a,- + 


1 


3 


49 


16440305941 


x 9 - 


2a; 8 - 


- 9a; 7 


+ 11a; 6 + 28a; 5 - 


18a; 4 - 


- 34a; 3 + 8a; 2 + 13a 


+ 1 


3 


229 


16898785417 


x 9 - 


-2a; 8 


-7a; 7 


+ 11a; 6 + 18a; 5 - 


17a; 4 


- 19a; 3 + 6a; 2 + 7a; 


+ 1 


1 


1 


16983563041 


X 9 - 


-a; 8 - 


- 8a; 7 


+ 7a; 6 + 21a; 5 - 15a; 4 - 


20a; 3 + 10a; 2 + 5a; 


- 1 


3 


361 


17515230173 


x 9 - 


-4a; 8 


-3a; 7 


+ 29a; 6 - 26a; 5 - 


24a; 4 


+ 34a; 3 - 2a,- 2 - 5a; 


+ 1 


3 


49 


18625670317 




x 9 - 


-9a; 7 - 


- a,- 6 + 23a; 5 + 4a; 4 


- 19a; 3 - 3a; 2 + 4a; + 1 




1 


1 


18756753353 




x 9 - 


3a; 8 - 


4a; 7 + 15a; 6 + 4a; E 


' - 22. 


r 4 - x 3 + 10a; 2 - 1 




1 


1 


19936446593 


x 9 


- 3a; ; 


s - 5a; 


7 + 17a; 6 + 7a; 5 - 


30a; 4 


- a,- 3 + 16a; 2 - 2a; - 


- 1 


3 


49 


20370652633 


x 9 


- 2a; ; 


s -8a; 


7 + 12a,- 6 + 15a; 5 - 


- 17x 


1 - 8a,- 3 + 8x 2 + x~ 


- 1 


1 


1 



References 

1. Julian Aguirre, Mikel Bilbao, and Juan Carlos Peral, The trace of totally positive 
algebraic integers, Math. Comp. 75 (2006), no. 253, 385-393. 

2. K. Belabas, A fast algorithm to compute cubic fields, Math. Comp. 66 (1997), no. 
219, 1213-1237. 

3. Manjul Bhargava, Gauss composition and generalizations, Algorithmic number 
theory (ANTS V, Sydney, 2002), Lecture Notes in Comput. Sci., vol. 2369, 
Springer, Berlin, 2002, 1-8. 

4. Henri Cohen, Advanced topics in computational number theory, Graduate Texts 
in Math., vol. 193, Springer- Verlag, New York, 2000. 

5. Henri Cohen and Francisco Diaz y Diaz, A polynomial reduction algorithm, Sem. 
Theor. Nombres Bordeaux (2) 3 (1991), no. 2, 351-360. 

6. Henri Cohen, Francisco Diaz y Diaz, and Michel Olivier, A table of totally com- 
plex number fields of small discriminants, Algorithmic number theory (ANTS 
III, Portland, Oregon, 1998), Lecture Notes in Comput. Sci., vol. 1423, Springer, 
Berlin, 1998, 381-391. 

7. Henri Cohen, Francisco Diaz y Diaz, and Michel Olivier, Constructing complete 
tables of quartic fields using Kummer theory, Math. Comp. 72 (2003), no. 242, 
941-951. 

8. Henry Cohn and Noam Elkies, New upper bounds on sphere packings I, Ann. 
Math. 157 (2003), 689-714. 

9. J.H. Conway and N.J. A. Sloane, Sphere packings, lattices and groups, 3rd. ed., 
Grund. der Math. Wissenschaften, vol. 290, Springer- Verlag, New York, 1999. 

10. Jesus De Loera, Raymond Hemmecke, Jeremiah Tauzer, and Ruriko Yoshia, Ef- 
fective lattice point counting in rational convex polytopes, J. Symbolic Comput. 
38 (2004), no. 4, 1273-1302. 

11. Jesus De Loera, LattE: Lattice point Enumeration, 2007, 
http : //www. math. ucdavis . edu/~latte/. 

12. Jordan S. Ellenberg and Akshay Venkatesh, The number of extensions of a number 
field with fixed degree and bounded discriminant, Ann. of Math. (2) 163 (2006), 
no. 2, 723-741. 



13. U. Fincke and M. Pohst, Improved methods for calculating vectors of short length 
in a lattice, including a complexity analysis, Math. Comp. 44 (1985), no. 170, 
463-471. 

14. Farshid Hajir and Christian Maire, Tamely ramified towers and discriminant 
bounds for number fields, Compositio Math. 128 (2001), no. 1, 35-53. 

15. Farshid Hajir and Christian Maire, Tamely ramified towers and discriminant 
bounds for number fields. II., J. Symbolic Comput. 33 (2002), no. 4, 415-423. 

16. Number field tables, ftp://megrez.math.u-bordeaux.fr/pub/numberfields/. 

17. KKiners-Malle, A database for number fields, 

http : //www. math. uni-duesseldorf . de/~klueners/minimum/minimum .html. 

18. Jiirgen Kliiners and Gunter Malle, A database for field extensions of the rationals, 
IMS J. Comput. Math. 4 (2001), 182-196. 

19. M. Kreuzer and H. Skarke, PALP: A Package for Analyzing Lattice Polytopes, 

2006, http : / /hep . itp . tuwien . ac . at/~kreuzer/CY/ CYpalp . html. 

20. A.K. Lenstra, H.W. Lenstra, L. Lovasz, Factoring polynomials with rational coef- 
ficients, Math. Ann. 261 (1982), no. 4, 515-534. 

21. D.D. Long, C. Maclachlan, A.W. Reid, Arithmetic Fuchsian groups of genus zero, 
Pure Appl. Math. Q. 2 (2006), no. 2, 569-599. 

22. Gunter Malle, The totally real primitive number fields of discriminant at most 10 9 , 
Algorithmic number theory, Lecture Notes in Comput. Sci., vol. 4076, Springer, 
Berlin, 2006, 114-123. 

23. Jason Martin, Improved bounds for discriminants of number fields, submitted. 

24. Jacques Martinet, Petits discriminants des corps de nombres, Number theory days 
(Exeter, 1980), London Math. Soc. Lecture Note Ser., vol. 56, Cambridge Univ. 
Press, Cambridge-New York, 1982, 151-193. 

25. Jacques Martinet, Tours de corps de classes et estimations de discriminants, In- 
vent. Math. 44 (1978), 65-73. 

26. Jacques Martinet, Methodes geometriques dans la recherche des petitis discrimi- 
nants, Sem. Theor. Nombres (1983-1984), Birkhauser-Boston, 147-179. 

27. A.M. Odlyzko, Bounds for discriminants and related estimates for class numbers, 
regulators and zeros of zeta functions: a survey of recent results, Sem. Theor. 
Nombres Bordeaux (2) 2 (1990), no. 1, 119-141. 

28. The PARI Group, PARI /CP (version 2.3.2), Bordeaux, 2006, 
http : //pari .math.u-bordeaux. f r/. 

29. Michael Pohst, On the computation of number fields of small discriminants in- 
cluding the minimum discriminants of sixth degree fields, J. Number Theory 14 
(1982), no. 1, 99-117. 

30. X.-F. Roblot, Totally real fields with small root discriminant, 
http : //math.univ-lyonl . f r/~roblot/tables .html. 

31. William Stein, SAGE Mathematics Software (version 2.8.12), The SAGE Group, 

2007, http://www.sagemath.org/. 

32. C.J. Smyth, The mean values of totally real algebraic integers, Math. Comp. 42 
(1984), no. 166, 663-681. 

33. Kisao Takeuchi, Totally real algebraic number fields of degree 9 with small dis- 
criminant, Saitama Math. J. 17 (1999), 63-85. 

34. Jan Verschelde, Algorithm 795: PHCpack: A general-purpose solver for polynomial 
systems by homotopy continuation, ACM Transactions on Mathematical Software 
25 (1999), no. 2, 251-276. 

35. John Voight, Totally real number fields, 

http: //www. cems .uvm.edu/~voight/nf-tables/. 



